3-CONNECTED CORES IN RANDOM PLANAR GRAPHS 



NIKOLAOS FOUNTOULAKIS AND KONSTANTINOS PANAGIOTOU 

Max-Planck-Institute for Informatics 
Saarbriicken, Germany 

Abstract. The study of the structural properties of large random planar graphs has become 
in recent years a field of intense research in computer science and discrete mathematics. 
Nowadays, a random planar graph is an important and challenging model for evaluating 
methods that are developed to study properties of random graphs from classes with structural 
side constraints. 

In this paper we focus on the structure of random biconnected planar graphs regarding 
the sizes of their 3-connected building blocks, which we call cores. In fact, we prove a general 
theorem regarding random biconnected graphs from various classes. If B„ is a graph drawn 
uniformly at random from a class B of labeled biconnected graphs, then we show that with 
probability 1 — o(l) as n — > oo, B„ belongs to exactly one of the following categories: 

(i) Either there is a unique giant core in B„, that is, there is a < c = c(B) < 1 such 
that the largest core contains ~ cn vertices, and every other core contains at most n a 
vertices, where < a = a(B) < 1; 

(ii) or all cores of B n contain 0(log n) vertices. 

Moreover, we find the critical condition that determines the category to which B n belongs, 
and also provide sharp concentration results for the counts of cores of all sizes between 1 
and n. As a corollary, we obtain that a random biconnected planar graph belongs to category 
(i), where in particular c = 0.765 . . . and a = 2/3. 



1. Introduction 

A fundamental discipline of computer science is the theoretical and practical evaluation of 
the performance of algorithms. From a theoretical viewpoint, many important computational 
problems turn out to be AA'P-hard or even very difficult to approximate, which leaves little 
hope for finding efficient algorithms that solve all instances. On the other hand, in practice 
it is widely observed that one can find satisfactory solutions efficiently, even when using 
algorithms that are known to perform very badly on certain inputs. In other words, it seems 
that "typical" instances are in some sense easier and the classical worst-case analysis may be 
too pessimistic. This motivates the study of the performance of algorithms from an average 
point of view. However, in order to perform such an analysis, it is necessary to specify an 
appropriate probability distribution on the set of all inputs. 

After having a suitable model at hand, an important ingredient of a meaningful average- 
case analysis is the precise knowledge of the typical structure of an input sampled from the 
specified distribution. For instance, in the context of graph problems, one wants to study the 
properties of graphs in the corresponding random graph model. However, this may be a very 
challenging task: if the space of inputs consists of graphs with global structural constraints, 
then the existence of dependencies can make the analysis formidable. 

A natural graph class with such constraints that has attracted the attention of researchers 
in computer science and discrete mathematics in recent years is the class of random planar 
graphs. In fact, since it was first studied by Denise, Vasconcellos, and Welsh [3], this model 
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has evolved as the primary example of studying random graphs from constrained classes, see 
e.g. PHIE]. From the perspective of the average-case analysis, it would be helpful to know that 
although a random planar graph inherits the dependencies that arise from the requirement of 
planarity, essentially most of its vertices behave independently of each other. However, the 
planarity condition makes almost all the tools and methods that have been used in the past 
decades for the analysis of classical random graphs models fail in the case of such a problem. 
Consequently, the development of new approaches is necessary and essential. 

One attempt to resolve this issue was taken recently by the second author and Steger 
The precise setting they considered is as follows. Let C be a class of labeled connected graphs, 
and let C n be a random graph from C with n vertices. The main idea in their work is to 
consider the maximal biconnected components of a connected graph C n . In this context, they 
showed, among other results, that under certain assumptions C n belongs a.a.s0 to exactly 
one of the following categories: 

(i) There is a unique giant biconnected component in C n . More precisely, there is a 
< c = c(C) < 1 such that the largest biconnected component contains ~ cn vertices, 
while every other component contains o(n) vertices. 

(ii) All biconnected components contain O(logn) vertices. 

Additionally, in [llj it was shown that random planar graphs belong to the former category, 
whereas e.g. random outerplanar graphs belong to the latter. Observe that for graphs that 
belong to category (ii), almost all pairs of vertices lie in different biconnected components, 
while this is not the case for graphs from the first category. A consequence of these facts is 
the following important observation. Random graphs from classes that belong to category (ii) 
"contain" in a well-defined sense plenty of independence. In particular, any such graph can 
be generated by choosing independently every one of its biconnected components, and gluing 
them together at the cut-vertices. As the biconnected components contain few vertices, and 
as they intersect each other only at single vertices, the impact of each block to the whole graph 
is small. So, such graphs resemble in a certain way the behavior of classical random graphs, 
where each edge is included independently with a specified probability, with the difference 
that here we choose the blocks independently of each other. However, random graphs from 
classes that belong to category (i) do not have this property: a lot of structure that we cannot 
control is "hidden" in the giant biconnected component, which contains a constant fraction 
of the vertices. 

This motivates a finer analysis regarding the typical structure of random biconnected 
graphs, which is the main topic of this work. In particular, we investigate how and un- 
der which conditions we can decompose a biconnected graph into building blocks of higher 
complexity. As we shall see shortly, such a decomposition is well-known and possible. We 
will show that we encounter again a fundamental dichotomy: depending on some critical 
condition, which we determine explicitly, we prove that large biconnected graphs have either 
only "small" building blocks of higher complexity, or a constant fraction of the vertices is 
contained in such a building block. Hence, we discover a picture that is completely analogous 
to the distribution of the sizes of the blocks in random connected graphs. 

Concerning the actual analysis, our methods not only differ, but also extend significantly 
some of the ideas presented in [11] with respect to two major issues. Firstly, essential ingre- 
dients of our proofs are precise asymptotic estimates for the number of biconnected graphs in 
the class in question. Here we resolve this issue in an analytic context and obtain very general 
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and precise enumeration results that are applicable in a wide variety of scenarios frequently 
encountered in modern theories of asymptotic enumeration. Secondly, regarding the combi- 
natorial aspect, the further decomposition of biconnected graphs in building blocks of higher 
complexity is considerably more involved than the decomposition of connected graphs into 
their biconnected components. So, although we use the same basic idea as in namely 
to sample algorithmically a graph from the class in question, the nature of the resulting al- 
gorithms is significantly more complex and cannot be analyzed with standard methods from 
probability theory. We resolve this issue by proposing a powerful method that relates by 
equation systems the simultaneous evolution of several random variables during the sampling 
process. We believe that this analysis opens up the possibility to study parameters of arbi- 
trarily complex classes that can be described within the framework of the symbolic method, 
which nowadays is a widely used tool to decompose various classes of combinatorial objects. 

Our results. The main idea of graph classification and enumeration that we will exploit 
goes back to Tutte |13j and consists of describing graphs in terms of building blocks of higher 
connectivity. For example, an arbitrary graph can be described by its connected components. 
Similarly, any connected graph decomposes uniquely into the set of its maximal biconnected 
components. For the class of biconnected graphs there is a similar decomposition with respect 
to 3-connected building blocks. Informally, one can give a recursive description of biconnected 
graphs in which 3-connected graphs have their edges replaced by biconnected graphs. This 
decomposition has been formalized through the notion of a network whose precise definition 
will be given in Section [2j The 3-connected building blocks which play the above role are 
called the cores of a biconnected graph. Using this decomposition, we can define classes 
of biconnected graphs by restricting the cores to be within a certain class of 3-connected 
graphs. For example, if we use the class of all 3-connected planar graphs as the base of 
the decomposition, then we obtain the classes consisting of all biconnected, connected, and 
general planar graphs. 

In this paper we study the distribution of the sizes of the cores in large random biconnected 
graphs, that are specified in Tutte's sense by their 3-connected building blocks. Our first 
result, in a very weakened form, is about random planar graphs. Here and in the remainder 
we write "x =L y" for the interval (x — y, x + y). 

Theorem 1.1. Let B n denote a random labeled biconnected planar graph with n vertices, and 
let u) = uo{n) be any function such that lim n ^ 00 u;(n) = oo. Let e > 0. Then, B n contains 
a.a.s. a unique largest core with (1 ± e)cn vertices, where c = 0.765... is given explicitely. 
Moreover, every other core contains at most n 2 l^uj(n) vertices, and there are (l±e)c^n cores 
with £ vertices, where 4 < £ = 0(n 2 / 3 ) and the ci are given explicitely. 

The main results of this paper are far more general. In fact, we show that a random planar 
graph is just a special case of a universal phenomenon. Let B n be a graph drawn uniformly at 
random from a class B of labeled biconnected graphs. In Theorems 12.31 and 12.41 we determine 
a critical condition according to which B n belongs a.a.s. to exactly one of those categories: 

(i) There is a unique giant core in B n : there are < c = c{B),a = a(B) < 1 such that 
the largest core contains (1 ± e)cn vertices, while every other core contains 0(n a ) 
vertices. 

(ii) All cores of B n contain O(logn) vertices. 

Moreover, in each of the above cases we determine for all 4 < £ < n the expected number 
of cores with £ vertices, and provide very sharp accompanying large deviation results. We 
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refer the reader to Section [2] for the formal presentation of our setting, and to Theorems 12.31 
and 12.41 for a precise formulation of the above statements. 

Our results have implications for several classes of biconnected graphs different from planar 
graphs. Let us mention selectively two examples. Writing Ex(Gi,... , G&) for the class of 
biconnected graphs that do not contain G\ , . . . , Gk as minors, the main result in [8] together 
with our result imply that a random graph from Ex(i^3 ] 3) on n vertices a.a.s. exhibits a 
behavior as above: it has a unique giant core whereas every other core is of lower order. Also, 
for the class ~Ex(K$ — e), which consists of all biconnected graphs whose cores belong to the 
set {-^3,3, K3 x K2} U {W n } n >3, where W n denotes a wheel onn + 1 vertices, again the main 
result in [8] together with our results imply that a random graph on n vertices from this class 
contains a.a.s. only cores having at most O(logn) vertices. 

Outline. As we mentioned above, the decomposition of a class of biconnected graphs is 
facilitated by the notion of a network, which we formally describe in the next section. This 
is a recursive decomposition that will allow us in Section [5] to design efficient algorithms 
that sample networks according to the so-called Boltzmann model. Such ideas were used 
for the first time e.g. in [I], where the number of vertices of a given degree in random 
dissections of convex polygons was studied. However, the nature of the decomposition in 
this work is very complex, and this makes the analysis of the samplers significantly more 
involved compared to the analysis in [T]. In particular, in order to sample networks, 
we use four different randomized algorithms, which call each other recursively in a nested 
and unpredictable fashion. This difficulty requires a completely new treatment as far as the 
analysis is concerned, and the details are presented in Section [5j 

Another crucial ingredient in our proof is a precise counting estimate for the number of 
networks, parametrized by the number of vertices and edges. To this end, we extend and 
further develop in Section certain known analytic tools, in order to make them applicable 
in our actual setting. Furthermore, they allow us to show a central limit theorem as well as 
tail bounds for the number of edges of a (general) random network with n vertices. Finally, 
in Section [7] we give the proofs of the main theorems and in Section [8] we give the proofs of a 
number of auxiliary results which are used in our arguments. 

Notational preliminaries. At this point let us introduce some necessary notation. Let C 

be a class of labeled graphs. For any positive integers n,m > 1, we denote by C n ^ m the 

subclass of C consisting of those graphs that have n labeled vertices and m edges, and we set 

C n = U m >oC njm . Moreover, we denote by C(x,y) the exponential generating function (egf) 

for C, where x marks the vertices and y marks the edges, i.e., the coefficient [x n y m ]C(x, y) 
I c l 

equals For a given y, we denote by pc(y) the singularity of C(x,y) with respect to x 

and, in particular, when y = 1 we write pc = Pc(l)- Also, we write Cb(y) := C(pc(y),y) 
and C := C(p c , 1). 

We now introduce a number of operations between combinatorial classes that are at the 
core of our analysis. For two graph classes X and y we denote by X x y the cartesian product 
of X and 3^ followed by a relabeling step, so as to guarantee that all labels are distinct. 
Note that the relation "A = X x y" expresses the fact that there is a bijection between the 
elements of A and pairs of elements from X and y, but it does not provide any information 
about how this bijection looks like, i.e., how to construct a graph in A from two graphs in X 
and y. The same is true for the operators described in the remainder. The class X + y 
consists of graphs that are either in X or in y. We denote by Set>k(X) the graph class 
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such that each object in it is an unordered collection of at least k graphs in X. Finally, the 
class X o e y consists of all graphs that are obtained from graphs from X, where each edge 
is replaced by a graph from y. Here we will usually assume that y is a class of graphs with 
two distinguished vertices, which simply means that we attach a graph G e € y at each edge 
e € E(G), where G € X, by identifying the special vertices of G e with the endpoints of e. 
This set of combinatorial operators (cartesian product, disjoint union, set, and substitution) 
appears frequently in modern theories of combinatorial enumeration, and it is beyond the 
scope of this work to survey them. For a very detailed description and numerous applications 
we refer to [6] and references therein. 

2. Networks 

Biconnected Graphs & Networks. Let B be a class of labeled biconnected graphs. In 
the remainder we assume that the graph consisting of a single edge is in B. Before we study 
the typical properties of a graph drawn uniformly at random from B n , that is, the subclass of 
B consisting of graphs on n vertices, let us introduce an auxiliary graph class that plays an 
important role in the decomposition of biconnected graphs. Following Trakhtenbrot [12] and 
Tutte |13j we define a network as a connected graph with two "special" vertices, called the left 
pole and the right pole, such that adding the edge between the poles the resulting (multi-)graph 
is biconnected. The remaining vertices of a network are called labeled vertices. Following 
standard notation, we will assume that in the egf of a class of networks the parameter x 
always marks the number of labeled vertices. 

The above description provides us with an explicit relation between the class B and the 
(corresponding) class of networks Af, which we shall describe now. Let B be the class contain- 
ing ordered pairs (B \ e, e) where B £ B and e G E{B), that is, an edge of B is distinguished 
and removed. Note that each B € B n , except from that where the underlying graph is a single 
distinguished edge, gives rise to two networks with n — 2 labeled vertices: one is obtained 
by removing the labels from the endpoints of the distinguished non-edge (and relabeling the 
remaining vertices such that only labels in {1, . . . , n — 2} appear), and the other is obtained 
by adding the distinguished non-edge to the underlying graph of B, removing the labels and 
performing a relabeling as before. Moreover, the single distinguished edge gives rise to the 
network that consists of two isolated poles. With X denoting the class consisting of a single 
labeled isolated vertex, and e being the network consisting of a single edge, this implies that 
the classes B and TV are due to the definition of Af related through 

(2.1) {Af + 1) x X 2 = (1 + e) x B. 

Note that this immediately translates to the relation 9B gy' v ^ = %- l+ i^f v ^ , which is the well- 
known dependency linking the egf's for B and (the associated) Af, see also e.g. |14] , 

Our Setting. The main setting in this work is the following. We will begin with a class of 
networks Af, which fulfills certain closure properties that are described below. This, combined 
with (|2.ip . defines a class of biconnected graphs B. Those B are precisely the ones that we 
will consider here. The following statement, along with our probability estimates, implies 
that the structural properties of random networks that occur asymptotically almost surely 
translate into almost sure properties of random biconnected graphs. Consequently, in the 
remainder of the paper we only need to consider random networks. Let V denote a graph 
property, that is, a set of graphs which is closed under automorphisms. With a slight abuse 
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Figure 1. Series and Parallel Networks. 

of terminology we say that a network N has a property V , if the graph that is obtained by 
putting labels to the poles has property V . The proof can be found in Section [8.1[ 

Proposition 2.1. Let M be a class of networks and let B be the class of the corresponding 
biconnected graphs. Moreover, let B n be a uniform random graph from B n , and N n a network 
that is drawn uniformly at random from M n . Suppose that P(N n _2 £P) > 1 — f(n — 2), 
where V is any property of graphs that is closed under automorphisms. Then P (B n € V) > 
1 — 2k f (n — 2), where nn is the maximum number of edges in a graph in B n . 

With the above facts at hand we are ready to define the classes of networks that we will 
consider. Let T denote a class of labeled 3-connected graphs that is closed under automor- 
phisms. Following [121 H3j. we define a class of networks with respect to T, denoted by M(T), 
inductively as follows (see also Figures [I] and [2]) : 

A network in M(T) is either an edge, whose endvertices are the poles, or is of type S (series 
network), or of type P (parallel network), or of type H (core network). 

Series networks: A network of type S consists of two networks N\ and N2, such that the 
right pole of Ni is identified with the left pole of N2. Here, Ni is restricted to be either an 
edge, or of type S or H, and N2 £ M(T). 

Parallel networks: A network of type P consists either of an edge and a non-empty set of 
networks, either of type S or of type H, where their right poles (left poles) are identified into 
a single right pole (left pole), or a set of networks of size at least 2, either of type S or of type 
H where the identification of the poles is as before. 

Core networks: Let T be the class of networks which are created by taking any graph in T, 
deleting an edge, and then turning its endvertices into poles. A network of type H consists 
of a network from T, where each edge is replaced by a network whose poles are identified in 
a unique way with the endvertices of the edges. 

In any of the above cases, we do the necessary relabeling in case there are conflicts, when two 
or more networks are joined through one of the above operations. With the above notation 
at hand, we can now introduce formally the classes of networks that we are interested in. 




Figure 2. A core network whose core is K4. 
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Definition 2.2. We say that M{T) is a-nice, if the following conditions hold. 

(A) For all y > 0, T(x, y) is analytic in C or has a unique dominant singularity at Pf{y) > 
0, and admits a local singular expansion of the form 



M n* >y) = E^)(i-^) 



where m £ N ; and the tk(y) are analytic functions. Furthermore, there exist 5, e > 
such that T(x,y) is analytic in a domain A = A(5, e) = {z : \z\ < pf(y) + 5, arg(z — 
Pf{y)) > e}. Moreover, the smallest integer k such that tfc(y) ^ and k/m N 
satisfies k/m = a. Also, the function pf{y) is strictly decreasing and continuously 
twice differentiable. 

(B) For y in a neighborhood of 1, let pn{u) be the radius of convergence of the egf enu- 
merating N{T) with respect to vertices and edges. Then 



The assumptions in the definition above are satisfied by several classes of graphs that 
have been studied in the literature, e.g. by series-parallel, planar or K^^-miaoi free graphs. 
However, we believe that the assumption (B) is redundant, in the sense that it follows from 
(A), but we are unfortunately unable to show it. The remaining conditions are minimal in the 
context of analytic combinatorics, as they form the "backbone" of essentially any counting 
result or limit theorem that can be derived. 

For notational convenience and in slight abuse of notation we will be writing T(x, y) for 
the egf T(x, y). With all the above definitions at hand we are now ready to present our main 
results. Given a network N, we denote by Ci(N) the number of vertices in a largest core 
of N. If £ is a positive integer and N is a network, we denote by c(£; N) the number of cores 
in iV which have £ vertices. More generally, if £ > 1, we denote by c(£, N) the number of 
cores in N whose number of vertices is at least £ and no more than ££. Set 



An important property of this function is that it defines implicitly the egf enumerating 
networks through the relation $>(x,y,N(x,y)) = 0; see Lemma 13.21 Our results imply 
that the sign of J^&(pn, 1, iVo) =: A dictates whether a random network has only at most 
logarithmically-sized cores (A > 0), or a giant core that contains a constant fraction of the 
vertices (A < 0). If u 6 {x,y,z}, we shall denote in the remainder by <& u (x,y, z) the partial 
derivative of &(x,y, z) with respect to u. Moreover, we denote by N n a random network 
from M n throughout and without further reference. In the case A > we show the following. 

Theorem 2.3. Let M{T) be an a-nice class for some a € M. \ {0}, and let 5, e > 0. Assume 
that & z (pn, 1j No) > and set r = j^fn^ ■ Then the following holds a.a.s. 

(i) For all 2 < £ < (1 — 5)log 1 / T n we have c(£; N n ) G (1 ± e)cin, where eg is given in 



Pjv(1) Piv(l) 



+ 



{pn(i)J ^ a 



(2.3) 
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Our next theorem deals with the case A < 0, where it is shown that the corresponding 
random networks have a.a.s. a unique giant core. 

Theorem 2.4. Let M{T) be an a-nice class for some a € M. \ {0}, and let 5, e > 0. Assume 
that & z (pn, 1,-Wo) < 0. Let u(n) be a function such that u(n) — > oo as n — > oo. Then the 
following holds a.a.s. 

(i) Let > be as in Lemma \7.4\ Then Ci(N n ) £ (1 ± e)7rn. 

(ii) For a// (nw(n)) 1/a < £ < Ci(N„) we feawe c(£; N„) = 0. 

/ \l/(a+l) 

(hi) For all 2 < £ < ( ^( n "i og n ) we ^ awe c (^5 N n ) S (1 ± s)cgn, where cp is given in 

Lemma \7.1\ 

(iv) Let e > 1. If2 <£ < i/ien c(*,#; N n ) G (1 ± e)c^n, w/iere is 

given in Lemma \7.6\ 

Moreover, for large £ the Equations (|7.5p . (|7.7p provide asymptotic expressions for cp andcg^g. 

Recently, Gimenez, Noy and Rue [8] obtained independently results regarding the asymp- 
totic distribution of the largest core in a random planar graph. Their methods are completely 
different and are based on techniques that are widely used in analytic combinatorics. However, 
in the present work, we obtain a fairly precise picture of the typical structure of a random 
biconnected graph whose 3-connected cores are taken from various general families, and we 
perform a very precise census of the smaller cores in all possible analytic regimes. 

Apart from the structural results concerning the sizes of the cores of N n , we also obtain a 
limit law for the number of edges of N n . 

Theorem 2.5. Let M{T) be an a-nice class of networks with respect to a class T . Let N n 
be a random network with n vertices. Then e(N n ) satisfies a central limit theorem, where 

the expected value is — °( n )> an< ^ ^ ae variance is proportional to n. Moreover, there 

exists £o > and a C > such that for all < e < Eq we have 

P(|e(N n ) - E(e(N n ))| > eE(e(N n ))) < exp (-Ce 2 n) . 
3. Preliminaries 

3.1. Boltzmann Sampling. The above decomposition of networks can be expressed in terms 
of certain combinatorial operators, which are omnipresent in modern theories of enumeration 
as well as random generation of combinatorial structures. We refer the reader to the book by 
Flajolet and Sedgewick [6] for a detailed exposition of these techniques. 

Duchon et al. [5] have shown that the above operations can be used to design Boltzmann 
samplers. If 7 is a graph, we denote by e(7) and v (7) the number of edges and the number 
of labeled vertices of 7, respectively. A Boltzmann sampler for Q is a randomized algorithm 
that generates any 7 € Q with probability G ^ x ^ x — • Note that if we set y = 1 and 
we condition on the order of the output graph being n, then the distribution of the output 
graphs is the uniform distribution on the set of graphs of order n. 

3.2. Networks. The decomposition of networks in Section [2] implies directly the following 
statement. 

Lemma 3.1. Let e denote the class of networks that consist of precisely one edge and X 
denote the class which consists of the empty graph with one labeled vertex. Moreover, let J\f 
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be an a-nice class of networks with 3-cores from the class T. Then the classes Af,S,V,Tt 
satisfy the following relations. 

N = e + S + V + H, 

S = ( e + V + H)xXxN, 

V = ex Set>i(5 + H) + Set> 2 (S + W), 

H = To e M. 

Using this fact we immediately obtain the following relation for the corresponding expo- 
nential generating functions. This result was shown, among others, in [121 114j . 

Lemma 3.2. The exponential generating functions enumerating networks satisfy 

N(x, y) = y + S(x,y) + P(x, y) + H(x, y), 

S(x, y) = {y + P(x, y) + H(x, y))xN(x, y), 

P(x, y) = (l + y) (e s ^+ H (^ - l) - S(x, y) - H(x, y), 
H(x,y) = T(x,N(x,y)). 
Moreover, N(x,y) satisfies the equation $>(x,y, N(x,y)) = 0, where $ is as in (|2.3p . 

4. Singularity Analysis for Networks 

A typical situation that is commonly encountered in modern scenarios of asymptotic enu- 
meration is the following (see also [B]) : we want to determine the precise asymptotic behavior 
of the coefficients of an unknown generating function F(x), which is given implicitly by some 
functional equation G(x,F(x)) = 0. Of course, only in very seldom cases it is possible to ob- 
tain F(x) explicitly, and generally, we have to resort to a completely different "program" to 
extract enough information from the given functional equation, so as to be able to understand 
the behavior of [x n ]F(x). 

The fundamental principle behind this program is that the knowledge of the behavior of an 
analytic function in the vicinity of its dominant singularities gives us very precise information 
about the asymptotics of its coefficients. So, we have to accomplish two tasks: first, to 
determine the dominant singularities of F(x), and second, to find the asymptotic expansion 
of our function near those singular points. The main objective of this section is to accomplish 
this program for the classes of a-nice networks. 

According to Lemma 13.21 the function N(x, y) is given by the solution of $(x, y, N(x, y)) = 0, 
where <3? is as in (|2.3p . Let y > be any fixed value. As N(x,y) has only non-negative 
coefficients, by applying Pringsheim's Theorem (see e.g. [6]), we infer that if the radius of 
convergence of N(x,y) is pn(v), then pN{y) is a singularity of N(x,y). 

Note that <3?(x,y,z) is not analytic at points that satisfy x > Pt(z), as T(x,z) ceases 
to be analytic there. Hence, we always have that Pn{v) < PT(No{y))- Moreover, suppose 
that for some pair (x',z') such that x' < Pt( z ') we have &(x',y, z 1 ) = and simultaneously 
<& z (x' ,y, z') 0. Then, by the Implicit Function Theorem, see e.g. [6], $ is locally invert- 
ible, which means that there exists an analytic continuation of N(x,y) around x = x' , and 
N(x',y) = z' . We can therefore obtain information about p by studying the distribution of 
the zeros of & z ( x i 2/> z )i i- e -; the points around which $ fails to be invertible. 

It turns out that the precise localization and nature of the singularity of N(x, y) is dictated 
by one of the following analytic conditions: 
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• Subcritical case: We have & z (pN(y), Vi No(y)) = 0. In this case, $ has a branch point 
at (pN{y),y, No(y)). Moreover, if $ was analytic at (p^f(y),y, iVo(y)), then it readily 
follows that Pn(v) < Pt(^o(u))- We further show in Lemma 14.21 that under certain 
mild assumptions the singularity of N(x, y) is of square-root type. 

• Supercritical case: We have & z (pN(y),y, No(y)) ^ 0. Here, the function <I> is not 
analytic at (pN(y),y, No(y)), and inevitably we have that pn{v) = Pt(^o(u))- In this 
setting, we determine the conditions (see Theorem 14.31 for a precise statement) under 
which the singularity type of some implicitly given function is "inherited" from the 
singularity of the implicit function. 

The following lemma will become very handy when we study the set of zeros of <£ 2 . 

Proposition 4.1. Let y > and suppose that a dominant singularity of N(x,y) is given 
by pn{v)- Let pe = pN(y)e ie . Then, for any < 6 < 2tt 

\<$> z {p e ,y,N(p e ,y))\ > -<$> z (p ,y,N(p ,y)). 

Proof. A simple calculation shows that & z (x,y,z) = T y (x,z) — ^q^n^^p • Let us abbrevi- 
ate N(pg,y) = Ng. By applying the inequality \a — b\ > \b\ — \a\, which is valid for all complex 
numbers a, b, we obtain 



\$z(po,y,N(p e ,y))\ > 



l-p g Ng(2 + peN ) 



\Ty(Pe,N(p e ,y))\. 



(l + N e )(l + Pe N d y\ 

As T(x, y) and N(x, y) have only non-negative coefficients, and T is aperiodic, the triangle 
inequality implies that 

\T y (pg,Ng)\ < T y (\pe\, \N(pg,y)\) < T y {\p e \,N(\p e \,y)) = T y (p ,N ). 

Similarly, we obtain that 

|1 - pgN e 2 {2 + pgNg)\ > 1 - |p fl iV fl 2 (2 + P gNg)\ > 1 - p N$(2 + po^o), 

and, 

\{1 + Ng)(l + pgNg) 2 \ < (I + N )(l + p N ) 2 . 

The proof completes by putting everything together. □ 

The remainder of this section deals separately with the subcritical and the critical case. 

The Subcritical Case. Our result in this section says that if $ ceases to be invertible at 
some analytic point inside its domain, then the egf enumerating networks has a singularity of 
the square-root type. 

Lemma 4.2. Let &(x,y,z) be as in (|2.3p . and letY be a compact subset o/(0, +oo). Suppose 
that for any yo £7 there exists a minimal Nq > and a minimal < xq < Pt(Nq) such that 

(4.1) $(x ,y ,N ) = and <S> z (x , yo, N ) = 0. 

Then there exists an e > such that N(x, y) admits a local representation 



N(x,y) = g(x,y) - h(x,y)Jl 



pn(vY 

for y £ Y, \x - p N (y)\ < e and \arg(x - pn(v))\ °. where g(x,y),h(x,y) and p N (y) are 

analytic around x = xq and y = yo, and Pn{vo) = xq and g(xo,yo) = Nq. Moreover, for 
any y £Y, N(x,y) is analytic in a domain A(p(y)). 
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Remark. Proofs of similar statements were found in the literature for e.g. the egf enumer- 
ating series-parallel networks [2]. However, the proof of the above claim requires more work: 
the function $ turns out to have in general negative coefficients, which makes the appli- 
cation of standard tools (see e.g. [1]) for this purpose impossible. Nevertheless, the proof 
that is presented here makes more or less use of standard techniques, and we include it for 
completeness. 

Proof. Let us first collect some basic properties of As T has non- negative coefficients, for 
any x, y > and z > 

_ , , m . , 1 + 7xz + 3x 2 z 2 + x 3 z 3 + 2x + 2xz 2 

$ zz {x, y, z) = T zz (x, z) + — — > 0, 

(1 + z) z {l + xz) 6 



and 



z 2 



(4.2) (a, y,z) = T x (x,z) + > 0. 

As Q(xo,yo, No) = & z (xo, Uo, No) = and & zz (xo, yo, Nq) / 0, the Weierstrass Preparation 
Theorem guarantees the existence of a function H(x, y, z), which is analytic around {xo,yo, No) 
and H(xq, yo, No) ^ 0, and two functions p(x, y),q(x, y) that are analytic around (xo, yo) and 
satisfy p(xo, yo) = q(xo, yo) = 0, such that 

(4.3) $(x, y, z) = H(x, y, z) {{z - N ) 2 + p(x, y){z - N ) + q(x, y)) . 

So, every function that satisfies $(x, y, N(x, y)) = is given in a neighborhood of (xo, yo) by 
one of the two functions 



JV(x, y) = A^o — ± y q{x, y). 

Note that q(xo,yo) = 0. Moreover, by differentiating both sides of (|4.3p by x we 

obtain that fe(^o,yo) / 0, as $ x {x , yo, N ) ^ 0. Hence Px( - X °^ q (xo,y ) / 0, and again 

the Preparation Theorem guarantees the existence of a function K(x,y), which is analytic 
around (xo,yo) and K(xo,yo) 7^ 0, and a function r(y), which is analytic around yo and 
r(yo) = 0, such that 

p(x, y) 2 

^ q(x, y) = K(x, y)(x - x + r(y)). 

By putting everything together we obtain that N has a local representation around (xo,yo) 
of the kind 

N{x,y) = g(x,y) ± h(x,y)Jl 



pn(vY 

where g,h are analytic around (xo,yo) and g{xo,yo) = No, and Pn(u) is analytic around yo 
and Pn(uo) = xo- Moreover, we choose h without loss of generality such that h(xo,yo) > 0. 

We will argue next that the "—"-sign above is the right choice. Denote the solution with 
the — sign by 2Vi, and the other one by N2. Note that by the minimality of xq and Nq we 
have $ z (x, y, N) ^ for < x < xo and < N < No such that <£(x, y, N) = 0, which implies 
that there is a unique analytic function that satisfies $(x, y, N{x, y)) = around 0. Moreover, 
note that the function N\ is increasing as x approaches xo from left, while N2 is decreasing. 
To show that the increasing branch is the function we are looking for, we show that there is 
a analytic curve that connects N\ to the function that is analytic around x = 0. 
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We show the claim by contradiction. Let y$ € Y. First, suppose that there is a < xq < x 
such that N(xo,yo) = Nq. Then $(xq, y, Nq) = $>(xq, yo, Nq) = 0, which is a contradiction, 
as <I> is monotone with respect to its first variable, due to (|4.2p . On the other hand, suppose 
that there is a Nq < iVo such that N(xo,yo) = Nq. Note that & z (xq, yo, 0) = — 1 < 0, and 
due to the minimality of Nq, for all < N < Nq we have that <& z (xq, yo, N) < 0. This again 
contradicts the assumption $>(xo, yo, Nq) = &(xq, yo, Nq). In other words, we have shown 
that the the solution of yo, N(x, yo)) = never "leaves" the rectangle {(x, N) | < x < 
xo,0 < N < Nq}, thus implying that that the increasing branch analytically continues the 
solution around 0. 

To complete the proof we will show that N(x, y) is analytic in a A-domain. By applying 
Proposition 14.11 we obtain for any pg = pN(y)e te , where < 9 < 2ir, that 

$ z (p e ,y,N(p e ,y)) > -$ z (p Q , y, N(p Q , y)) = 0. 

Thus, the Implicit Function Theorem shows that there is an analytic continuation of N(x, y) 
around p$, where < 9 < 2ir. This completes the proof. □ 



The Critical Case & Transfer of Singular Expansions. The following result provides 
us with a local expansion of a function / that is given implicitly by a functional equation 
/ = G{x,y,f). In contrast to the previous result, here we handle general functions G at 
points where they are not analytic, thus extending a result from [4], where just the case of 
singular behavior of type 3/2 was treated. This result will be used in the proof of Lemma l4~4"l 
which complements the statement of Lemma [4. 2 1 in the case where the function $ is invertible 
in its whole disc of convergence. 

Theorem 4.3. Suppose that G{x, y, f) has a local representation of the form 

j ^ k/m 



G(x, y, f) = g(x, y, f) + h(x, y, f) 1 



r{x,y) 



where k,m 6 N, 1 < k/m £ N, the functions g,h,r are analytic around (xQ,yo, /o), and 
satisfy 

9f(x ,yo,fo)j^l, h(xQ,yQ,fo)^Q, r(xQ,y )^0, r x (x , y ) ^ g x (x , y , f ). 

Suppose that f = f(x,y) is a solution of the functional equation f = G(x,y,f) such that 
f{xQ, yo) = /o- Then f admits a local representation of the form 

k/m 

(4.4) f(x,y)=g(x,y)+h(x,y)[l 



p{y), 

where g,h,p are analytic at (xo,yo) and h(xo,yo) / and p(yo) = xq. 

The next lemma takes care of the cases in which $ is not analytic at (pN(y),y, N(pN(y),y)). 

Lemma 4.4. Let &(x,y,z) be as in (|2.3p . and letY be a compact subset o/(0,+oo). Suppose 
that for all yo € Y and all Nq > and < xo < Pt(Nq) such that <&(xo,yo, Nq) = it holds 

<S> z {xo,yo,No)^0. 
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Then, ifT{x,y) admits a local representation of the form ()2.2p and is analytic in a A-domain, 
there exists an e > such that N(x, y) admits a local representation 

k/m 

N(x,y) = g(x,y) + h(x,y) ( 1 



Pn{v) 

for y E Y, \x — Pn(v)\ < £ and \arg(x - Pn{v))\ / 0, where g{x,y),h(x,y) and p N (y) are 
analytic around x = xq and y = yo, and g(pN(yo),Uo) = Nq. Moreover, for any y EY, N(x, y) 
is analytic in a domain A{p^{y)). 

Proof. First of all, note that the assumption on T(x, y) implies that there are two functions 
g(x,y,z) and h(x,y,z), which are analytic at (xo, yo, No), such that locally 

/ x \ k/m 
$(x,y,z) = g(x,y,z) + h(x,y,z) 1 -- 

V PT(Z)J 

Our aim is to apply Theorem 14. 31 in order to determine a local expansion for N(x,y), which 
is implicitly defined by §{x,y,N{x,y)) = 0. To achieve this, we shall first show that we can 
"switch" in the above expression between local expansions in terms of x and z. Moreover, we 
further will check the remaining conditions in Theorem 14.31 

Recall that p? is a strictly decreasing function that attains only positive values. Hence we 
have that pt(No) ^ and p' T (No) ^ 0. The Weierstrass preparation theorem implies that 
there is a function H{x, z), which is analytic at (xo, Nq) and H{xq, Nq) ^ 0, such that 

Pt{z) — x = H(x, z)(r(x) — z), 

where r{x) is the analytic inverse function of px in a neighborhood of xq. Hence we infer that 
there is an analytic function h{x, y, z) such that 



V, z) = 9(x, y, z) + h(x, y,z) [l 



, k/m 
Z x ' 



r(x] 



It is now routine to check from our assumptions the remaining preconditions of Theorem 14.31 
Hence, we obtain the claimed local representation for N. 

To complete the proof we show that N is analytic in a A-domain. For the sake of con- 
tradiction, suppose that there is another singularity on the circle of convergence of N(x,y), 
and write pe = pN(y)e t0 , where < 6 < 2ir. Note that $ is analytic at (p(y),y, N(p(y),y)), 
because T is analytic there. So, we must have & z (p(y),y, N(p(y),y)) = 0, as otherwise the 
Implicit Function Theorem would guarantee the existence of an analytic continuation of N 
around p$. Then, by applying Proposition 14.11 we obtain 

(4.5) = $ z (p e , y, N( Pe ,y)) > -$ z {p ,y, N( Po ,y)) 

Note that c &(0, y, y) = 0, and hence N(0,y) = y (this follows also from the definition of 
networks, as there is precisely one network that consists of a single edge). A straightforward 
calculation shows that & z (0,y, y) = —j^j < 0. Hence, due to our assumptions we have for 
all < x < xq that & z (x, y, N(x, y)) < 0, which implies that (|4.5j) is a contradiction. □ 
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Proof of Theorem \2.5[ Let y be in some fixed complex neighborhood of 1. By combining 
Lemma [4.21 with Lemma [4.41 we obtain that there are analytic functions g, h such that locally 
around Pn(u) 

(4.6) N{x,y) = g(x,y)+h(x,y) (l 



pn{v), 

where r = 1/2 if <&(p N (y),y, N {y)) > 0, and r = a if $(p N (y),y, N (y)) < 0. By applying 
the Transfer Theorem (see Corollary VI. 1 in [B]) we obtain that the probability generating 
function p n (u) for the number of edges in N n satisfies 

p n (u) = h(x,u)——p N (u)- n (l + 0(l)). 

r(r) 

The proof finishes by applying the Quasi-Powers Theorem (see Theorem IX. 8 in [6]), and by 
exploiting Propoerty (B) in the definition of nice classes of networks. □ 



5. Sampling Networks 

Although the existence of Boltzmann samplers for several classes of networks is guaranteed 
by the work of Duchon et al. [5], we shall nevertheless describe them explicitly in this section. 
This will not only allow us to introduce some necessary notation that will be used in the rest 
of the paper, but it will also make possible to relate the number of 3-cores in networks to 
certain random decisions performed by the samplers. We will make such statements more 
precise in Section 15.21 

5.1. Boltzmann Samplers for Networks. Using the rules that are described in Section ETT1 
we construct the Boltzmann samplers for several classes of networks. 

We begin with the class TV. Note that any network in J\f is either an edge, or an 5-network, 
or a P-network, or a TL network. Hence, a Boltzmann sampler for J\f will call a sampler for a 
subclass with probability proportional to the value of the generating function of this subclass. 
More precisely, we say that a variable X is network- distributed with parameters x and y, X ~ 
Net(x,y), if its domain is the set of symbols Net = {e, S, P, H} and for any s € f2 Net it 
holds ¥(X = s) = jyf'^l . Then the sampler for J\f can be described concisely as follows. 

TN(x,y) : s <- Net(x,y) 
return Ts(x, y) 

Next we describe the sampler for S. For two networks N\ and N2 we will write N1N2 
for the network that is obtained by identifying the right pole of N\ with the left pole of A^2- 
Recall that S = (e+V+Tt) xM. Hence, FS(x, y) has to choose among the classes e, V, and H 
with the right probability. More precisely, we say that a variable X is series-distributed with 
parameters x and y, X ~ Ser(x,y), if its domain is the set of symbols f2 Ser = {e, P, H} and 
for any s 6 r2 Ser it holds F(X = s) = g^'^ . Then the sampler for S is given by the following 
procedure. 

TS(x,y) : s <- Ser(x,y) 
£ <- s(x,y) 
r <- N(x,y) 

return £r, relabeling randomly its vertices 



3-CONNECTED CORES IN RANDOM PLANAR GRAPHS 



15 



In the sequel we describe the sampler for V . A parallel network either consists of an edge 
and a set of at least one S or Ti network (= Set>i(<S + Ti)), or it consists of a set of at least 
two iS or Ti networks (= Set>2(<S + TC)). This implies that TP(x, y) has to first to choose one 
of the two possibilities with the right probability, and then to sample a set (with a given lower 
bound on the number of elements) from S + Ti according to the Boltzmann distribution. 

Let us introduce some notation before we describe formally the sampler. We say that 
a variable X is parallel- distributed with parameters x and y, and write X ~ Par(x,y), if 

X ~ 1 + Be( e — V) p{x y)^'^ H ^ x ' yS> ). Finally, we say that a variable is sh- distributed 
with parameters x and y, X ~ sh(x,y), if its domain is the set of symbols J7 sh = {S, H} and 
for s G f2 S h ^ holds ¥(X = s) = s ^ x ~rqrfj7^ y ^ ■ Then TP works as follows. 

TP(x,y) : p^Par(x,y) 

k <- Po> p (S(x, y) + H(x, y)) 

for i = 1 . . . k 

bi <- sh(x,y) 
Pi <- Tbi(x,y) 

construct a network P by identifying the left and right poles of p\, . . . ,pi- 
relabel randomly the vertices of P 

if p = 1 then return P, where the poles are joined by an edge 
else return P 

Finally, we describe the sampler for Ti. Recall that a W-network is obtained by substituting 
the edges of some graph from T by graphs from J\f. Here we will assume that we have 
an auxiliary sampler TT(x,y), which samples graphs from T according to the Boltzmann 
distribution. Then the sampler for Ti can be described as follows. 

TH(x,y) : T <— TT{x, N(x, y)) 
foreach edge e of T 
le ^YN{ Xl y) 

replace every e in T by 7<= 

return T, relabeling randomly its vertices 

This completes the definitions of all samplers that we shall exploit. From the theory 
developed in Duchon et al. [5] we immediately obtain the following statement. 

Lemma 5.1. Let N € M such that N has n labeled vertices and m edges. Then ¥(TN(x, y) = 
T\N(x,y) ■ The same statement is also true for networks in the classes S,V,Ti and the 
corresponding samplers given above. 

Using the definition of networks we obtain readily the following property. 

Lemma 5.2. Let M{T) be an a-nice class of networks such that a£l \ {0}. Let 

/ 5/2 , tf$ g (p N ,l,N o )>0 
P \ a , if&,(p N ,l,N o )<0 ' 

Then there is a b > such that F(TN(p N , 1) G M n ) ~ bn" 13 . 



Proof. The proof follows from 14.61 for the special case y = 1 and the Transfer Theorem 
(Corollary VI. 1 in [6]). □ 
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5.2. Combinatorics of the Samplers. Although the Boltzmann samplers are randomized 
algorithms, we are going to analyze them as if they were deterministic algorithms. In par- 
ticular, we are going to assume that all their parameters are read from lists, which contain 
independent random samples of the appropriate parameters. 

Let ^et be an infinite list of values from n Net . Similarly, £ser> ^Par and ^ S h denote infinite 
lists of values from Ser , f2 Par and O sh , respectively, and £j list of integers > j for j € {1,2}. 
Finally, let £q- be a sequence of graphs from T. 

Note that if we choose the values in the lists according to the corresponding probability 
distributions (i.e., the values in £ Net according to independent network-distributed variables 
with parameters x and y, the values in ^ser according to independent series-distributed vari- 
ables, . . . ), then the probability that the deterministic counterparts of the samplers generate 
a network is equal to the probability in the Boltzmann model. 

Using this observation, the proofs of our main results proceed according to the following 
schema. 

(1) Relate properties of a network generated by a sampler to properties of the values in 
the lists 4et,4er,4ar, 4hj 4, h and It- 

(2) Show that the desired properties are observed with high probability, if the values in 
the lists are chosen independently according to the corresponding probability distri- 
butions. 

The next two statements take care of step (1). More precisely, the following statement gives 
us the number of 3-cores in a network generated by TN. Its proof is straightforward and 
therefore omitted. 

Lemma 5.3. Suppose that the network N was generated by TN by using the first aj- graphs 
in It ■ Then 

c(k;N) = \{l < i < a r \ v(£ r [i]) = k -2}| . 

The next claim gives us some general relations among the values used by the sampler, which 
will be very useful later. The proof can be found in Section 15.21 and can be performed by 
looking closer at how the sampler constructs a large network out of smaller ones. We denote 
by 1(E) the indicator variable for the event E, i.e., 1(E) = 1 if E occurs, and 1(E) = 
otherwise. 

Lemma 5.4. Suppose that the network N was generated by TN by using the first A s values 
in l s , where s £ {Net, Ser, Par, sh, 1, 2, T}. Then the following statements are true. 

(5.1) v(N) = A Sei + <^[i]), 

l<i<A T 

(5.2) e(N) = Y, 1 (4etW=e)+ £ l(£ Sei [i\ = e) + £ l(£ Pax [i\ = 1), 

l<i<A„ et l<i<A Ser l<i<A P ar 

(5.3) A Net = l + A Ser + Yl < l r[i]), 

l<i<A T 

(5.4) Aj = J2 l(Wi]=j), 

l<i<A Par 

(5.5) A Par = Yl l(4etM=P)+ Y 1(W*1 = -P)- 

l<i<A Net l<i<A S er 
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Moreover, 

(5.6) A sh 

(5.7) A Ser 

(5.8) A T 



E 4W+ E « 

l<i<Ai l<i<A 2 

E l(4«t[»] = 5)+ E l(4hW = S), 



Ki<A« 



Ki<A s 



E lM1=fl)+ E 1(W*J=#) + E 1(M*]=#)- 



l<l<A Ne t l<i<A S er l<j<A sh 

6. Cores In Random Networks 

In our proofs we will be using repeatedly the Chernoff bound on the probability that a 
binomially distributed random variable deviates significantly from its expected value. We 
will use the version which appears in [9]. In particular, let X be a binomially distributed 
random variable and let t> 0. Then 

/ t 2 

(6.1) P(|X-E(X)| > t) < 2exp — — - 

v ; Vl v ;| '- v \ 2(E(X)+i/3), 

Another technical ingredient in our proof is the following lemma. It gives a Chernoff-type 
bound for the sum of independent Po>j(/x) variables, and its proof uses exponential generating 
functions. For the sake of completeness, it can be found in Section [8.31 

Lemma 6.1. Let X\, . . . , X r be independent Po>j(p) variables, where p > and j € {1, 2}. 

Let Y r i = Xa<i<r' ^-i, where < r' < r. Then, there is aC > such that for any < e < 1 
and sufficiently large r 

P (3 < r' < r : \Y r > — E (Y r /)| > er) < e _Ce2r . 

In the following we denote by A s the random variable counting the number of values 
from £ s used by an execution TN(pN(y),y), where s € {Net, Ser, Par, sh, 1, 2, T}. More- 
over, we write Vr and Eq- for the total number of labeled vertices and edges in all cores 
of TN(p N (y),y), i.e., 

At A t 
V T = J2v(£ T [i]) = ^2(k-2)c(k;TN(p N (y),y)), and E T = Y^e{i T [r]). 
i=l k i=l 

Finally, let us denote by N the output of the Boltzmann sampler TN(p^, 1). 

Lemma 6.2. Let < e < 1. There is a constant C > such that for sufficiently large n and 
any Z G {Aj et , A 5ei , A PaT , V r ,E T } 

P (|Z - zn\ < en | N e A/" n ) > 1 - e" Ce2n , 

and 2 € {aNeti osen «Par, ^Ti er}, where a = [a Net , a Ser , a Par , t>r, er] T unique solution 

of the linear system Ma = r, with 



(6.2) 
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Proof. First, by applying Lemma 15.41 Statements (|5,ip and (|5.3p . we obtain that for every 
point in the conditional probability space "N E Af n " 

n = A SeT + Vj and A Net = A SeT + Er, 

from which we immediately obtain that 1 = + — and = + — . These two facts 
are the second and the last line of the linear system above, and thus hold with probability 1. 
In the remainder we argue that all other equations are true with probability at least 1— e~ Ce n . 
This completes the proof of the lemma with the following reasoning. The determinant of M 
equals 

D = PN N$ + (p N + 1)N + 1. 

Now, since N(x) has only non-negative coefficients, we infer that D / and the proof is 
finished. 

Let us continue with an auxiliary observation. Consider e.g. Statement (|5.4p in Lemma [5 .41 
Our aim is to translate this statement into a high probability statement for the relation of the 
random variables Aj and A PaT . For this, let pj = P(Par(pjv,l) = j) and note by applying 
Lemma 15.21 we infer that there is a c > such that for large n 

P (Aj (1 ± e)p j A Pax ± en | N £ J\f n ) < cn^P (Aj g" (1 ± e)pjA Pax ± en) . 

By applying (|5.4p we obtain Aj = Ylt={ l(^Par[*] = j)- Using this, we infer that the proba- 
bility of the above event is at most 



cn 



(bL > 1 : l(4ar[»] = j) t (1 ± e) Pj L ± en^J 



and a simple union bound together with the Chernoff bounds imply that there is a C > 
such that 

P (Aj g" (1 ± e) Pj A Par ± en | N £ M n ) < e~ C£2n . 
In other words, we have demonstrated that ()5.4[) translates into the statement 

(6.3) Aj £ (1 ± e)pjA Par ± en 

with probability > 1 — e - ^ 2 " - . Now, with exactly the same line of reasoning we can infer 
from (|5.2p and (|5.5p that if we condition on "N £ A/" n ", with probability at least 1 — e~ Ce n 

( 1 o\rNn e S +H _ ^ \ 

(6.4) e(N) £ (l±e)(—A Net + ^-^A Ser + A Pai ± en, 

(6.5) A Par £ (l±e)(^A Net + ^^A Se ^±en. 

Moreover, by applying Lemma 16.11 instead of the Chernoff bounds we infer from (|5,6p that 
with probability at least 1 — e~ Ce n 

rfifi N a - n , , f (So + gp)e g ° +g ° Qgg + ff )(e g o+^ -1) \ 
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Finally, from (|5.7p and (|5.8p we infer again by the Chernoff bounds that with probability at 
least 1 - e~ Ce2n 

(6.7) A Ser G (l±e) f^Net+ Q S ° U A sh ) ± en, 

fa q\ A r- /i _l \ f H ° A , PnHoNo H \ 

(6.8) A T G (l±e) TT^Net H 7j ^Ser + TT~ TT^sh ± 

\iVo Oo Oo + -Ho / 

Moreover, Theorem 12.51 implies that there is a B > such that 

P (e(N) G (1 ± e)/m | N G An) > 1 — e~ Be2n . 

Let X = ^[Ajet> ^4ser, ^.pan Vt, Eq-] T ■ By combining all the above facts together, and using 
the fact N(x, y) = y + (1 + y)(e s ( x ' y ^ +H ( x ' y } — 1), which follows from Lemma l3T2l we infer that 
with probability at least 1 — e~ c £ n there is a c > such that 

(1 +e)M ■ X > r - e ■ c and r + e ■ c > (1 - s)M ■ X 

where c = [c, c, c, c, c] T . The proof completes by elementary linear algebra algebra, and by 
choosing the e in this proof as e/c', for a suitable d > 0. □ 

Note that the above statement does not yield any information about the total number of 
cores in a random network with n vertices. However, with a little additional work we arrive 
at the following result. 

Corollary 6.3. Let < e < 1. Then, there is a C > such that 

P[\A T - a r n\ > earn \ N G jV n ] < e ~ Ce2arn , 
where a r = 2fiT(p N , N ) and fi = -^tjj ■ 

Proof. By solving (|6.2[) we obtain explicit (but lengthy) expressions for the highly probable 
values of Aj et , ^4ser> and A sh . Then, by using (|6.8p we obtain after elementary algebraic 
manipulations the claimed statement. □ 



7. Proof of Theorems 12.31 and 12.41 

In this section we perform the proof of our main results. We shall denote throughout by N n 
a network drawn uniformly at random from M n and by N a graph generated by TN(p^, 1). 

7.1. Small Cores. As a first application of Corollary 16. 31 we prove that the counts of "small" 
cores in N n are sharply concentrated around a specific value. Before we proceed let us make 
an auxiliary observation that will be used several times. Let M{T) be an a-nice class of 
networks. Let B n C J\f n be some (bad) property of networks, and let N G M n . Moreover, let 
Aq- = Aq-(N) be the total number of cores in N. As the random network N has the same 
chance of being any N G M n we obtain by applying Corollary 16.31 that there is a C > such 
that 

P(N n G B n ) = P(N G B n I N G M n ) < P(N G B n , A r € (1 ± e/2)a T n \ N G An) + e~ c ' e2n , 
where ar = — 2 ^|^ T(pjv, Nq). Applying Lemma I5T2"1 we deduce that 
(7.1) P(N n G B n ) < 0(nP) ■ P(N G B n , A r G (1 ± e/2)a r n) + e"^ 2 ™. 
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The parameter (3 is as in Lemma l5.2i In particular, j3 = 5/2 if «£ z (p./v> 1, Nq) > 0, and /3 = a 
otherwise. This is a fact that we shall use below several times. We start by counting cores of 
a given size. 

Lemma 7.1. LetJ\f(T) be an a-nice class of networks for some a G M\{0} ; and letO < e < 1. 
Let aq- be defined as in Corollary 1 6. 31 (3 as in Lemma \5.Sh and define for k > 3 the quantities 

(7.2) p k := 2] J. (X '^ ) :^ 2 and k : = max{£ | p e n > VI e~ 2 a? {(3 + 1) log n}. 
T(p N ,N ) 

Then there exists an constant C = C (J\f) > such that for large n and k < ko 

¥(c(k; N„) G (1 ± e) a r p k n) > 1 - e' ^. 

Proof. Let B n C J\f n be the set of networks whose number of cores with precisely k vertices 
is not in (1 ± e) arPk^- By applying (|7.ip we infer that it sufficient to show that, say, P(N G 
B n , A r G (1 ± e/2)a r n) < e - c '^P^ n . 

Note that Lemma [531 asserts that c(fe; N) = Ylt=i l(f(^TH) = k — 2). As every T G T has 
vertices iff it has k — 2 labeled vertices, we infer also that p k = P(u(#r[*]) = k — 2). So, 

P (N G B R , A r G (l ± |) fl T n) < P (bL G (l ± |) a r n : l(v(£ T [i]) =k-2)#(l± e)a r p k n 

1 2 

By the Chernoff bounds and a union bound this probability is, say, at most 2ne~T6 e aT P kn . 
The proof then completes with (|7.1|) and the choice of ko for large n. □ 



With the above Lemma at hand we are ready to prove the first statement in Theorem 
and the third statement in Theorem 12.41 Note that the analytic properties of T asserted in 
Definition 12.21 imply with the Transfer Theorem (see e.g. Corollary VI. 1 in [6]) that 

(7.3) [x k - 2 ]T(x, No) ~ t^iMk— V(iV )- fe+2 . 

Proof of Theorem \2.3\, (i). The condition ^ z (p^,l, Nq) > together with Lemma 14.21 imply 
that pn < Pt(No). Moreover, the definition of p k in (17.2j) and (I7,3j) assert that there is a 
C > such that 

(7.4) p k ~ Ck~ a - X ( P ^ \ = Ck- a ^r k . 



p T {N 

As r < 1, we infer that we can apply Lemma 17.11 for ko = (1 — 5) log 1 ^ T n whenever n is 
sufficiently large. This completes the proof. □ 



Proof of Theorem \2.4\ (Hi). Here, the condition Q z {pn-, lj-^o) < together with Lemma 
imply that pn = Pt(Nq). The definition of pk in (I7.2p and (|7.3f) assert that there is a C > 
such that 

(7.5) p k ~ Ck^ 1 (— =Ck- a -\ 



PT (N ) 

We infer that we can apply Lemma 17.11 for ko = ( ^"logrJ 1 ^"^ wnenever n is sufficiently 
large. □ 
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The next lemma deals with cores that contain more than ( , 1\ vertices. The 

v uiyn) log n > 

proof is essentially the same as the proof of Lemma 17. II and hence omitted. 

Lemma 7.2. LetJ\f(T) be an a-nice class of networks for some a € M\{0} ; and letO < e < 1. 
Let aq- be defined as in Corollary \6. 31 (3 as in Lemma \5.2[ and define for k > 3 and £ > 1 

(7-6) p k ^ k : = V] — Ln/^' °\ Pn and k := max{£ \ pi^n > 17e~ 2 a 7 ^ 1 (/3 + 1) log n}. 

f^ k T{p N ,N ) 

Then there exists an constant C = C (J\f) > such that for large n and k <k$ 

F(c(k,Zk; N„) £ (lie) a rP k^n) > 1 - e~ c ^^ n . 

Proof of Theorem \2.4\ (in). The condition & z {pn-, 1 5 Nq) < together with Lemma 14.41 imply 
that pn = pt(Nq). The definition of Pk£k in (17.6j) and (IT. 3f) assert that for £ > 1 there is a 
C > such that 

(7.7) M ~ Cl ~ a ~ 1 ~ ^""C 1 " r*)- 

We infer that we can apply Lemma [7.21 for ko = ( ^( n )" ogn ) 1//ct whenever n is sufficiently large. 

□ 

We now consider the case Q z {pNi Ij-^o) > 0- The next statement deals with the cases in 
Theorem 12.31 not covered by Lemma 17. 11 

Lemma 7.3. Let N{T) be an a-nice class of networks for some a£l \ {0} ; and let e > 0. 
Assume that ® z (p N (y),y, N (y)) > 0, and let r = p ^ N(>) ■ Then 

P (Ci(N n ) > (5/2 + e) log 1/T nj = (n~ e ). 

Moreover, we have 

c ^(1 - e) log 1/r n, ^ log 1/T n; N n ^ > n 2 ^ = o(l). 

Proof. Let 23 n C A/" n be the set of networks in which the largest core has size > (5/2 + 
e) logi/ T n. By applying (|7.ip we obtain 

P(N n € B n ) < 0(n 5 / 2 ) • P(N G £„, Ar6(l± e/2)a r n) + e" c ' £2n . 
Note that Lemma 15.31 asserts that 

NGB n => £l(w(4r[i]) > (5/2 + e)log 1/T n-2) > 0. 

i=l 

As P(u(£r[i]) = fe) = — ^feS^ 2 w e obtain by using {73]) that there is a C > such 
that for large n 

v(£ T [i]) > (5/2 + e) log 1/r n - 2) < C ^ k~ a ~ l T k . 

fc>(5/2+ £ )log 1/T n-2 
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The condition $ z (pN,l, N ) > together with Lemma [4.21 imply that pjy < pj>(N ), and 
hence r < 1. Thus, P (v (#r[*]) > (5/2 + e) log 1 / T nj = o(n~ 5 / 2_£ ). The proof of the first 
statement completes with Markov's inequality. 

The second statement follows similarly by estimating the probability 

;i - e) log 1/T n < v(£t[{\) < (5/2 + e) log 1/r n) 

with the same technique as above. The straightforward details are left to the reader. □ 

7.2. The Largest Core. In this section we prove the first two statements in Theorem 12.41 
In particular, we show the following. 

Lemma 7.4. Let M{T) be an a-nice class of networks for some a£l \ {0}, and let e > 0. 
Assume that & z (pi\[(y),y, No(y)) < 0, and let u(n) be a function such that linin^oo uj{n) = oo. 
Then, a.a.s. 

|Ci(N n ) - 7rn| < en, 

where 77- := vq- — a T^rfp^I^f^ ' an< ^ V T ^ s as * n Lemma \6.2\ and aq- as in Corollary 1 6. 31 
Moreover, there are a.a.s. no other cores in N„ with more than u(n) ■ n 1//a vertices. 

Proof. Set no := uj(n) ■ n L l a . We will use a counting argument to show that the number of 
networks that have a unique core with more than no vertices is asymptotically much larger 
than the number A2 of networks that have at least two such cores. 

The singularity expansion of N(x, y) with respect to x (Lemma 14. 4p together with the 
Transfer Theorem (Corollary VI. 1 in [6]) imply that there exists a constant c such that 
as i — ► 00 

(7.8) \NA = {l + o{l))cp- N H- a - 1 t\. 

We bound A2 as follows. Let TV be a network having n labeled vertices. Assume that iV 
contains at least two cores each having at least no vertices. Then there is either a cut-edge 
that splits N in two networks that contain at least no vertices each, or such a splitting can be 
obtained if we choose the two poles as the cut-set. So, every such network can be described 
by a triple (Ni,e, N2), where e G N\ (or, in slight abuse of notation, e contains the poles 
of Ni), v(N\),v(N2) > no — 2, and we can construct N by identifying the poles of ./V2 with 
the endpoints of e. 

Note that by Theorem 12.51 we can assume that with probability 1 — o(l) there are at 
most Kn choices for e, for some constant K > 0. Moreover, there are (^/^-j) ways to choose 
the labels for N\. By summing over all choices for v(Ni) we obtain for n large enough that 



A 2 < 2Kn ]T r j |A4| \X n -s\ + o(\Af n \) 



no/2<s<n— no 

< 4ifn Yl (") • PN s ~ a ~ ls] ■ PN n+S (™ ~ sY a -\n - s)\ + o(\M n \) 

n /2<s<n-n 

= AKnp- N n n\ s~ a - l {n- S y a - l +o(\N n \) 
no /2<s<n—no 

= SKnp- n n\ Y s" Q " 1 («-s)- Q - 1 + o(K|). 

no/2<s<n/2 
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The last sum can be bounded above as follows: 

n /2<s<n/2 n /2<s<n/2 " 

= O (n _a_1 no a ) = O (n" -^^)- 1 ) . 

So using (|7.8p with £ = n, we obtain = o(|jV n |). 

This shows that there is a unique largest core with more than uq vertices. To complete 
the proof we show the claimed bound on its size. First, by applying Lemma 16.21 we infer that 
a.a.s. the total number ^fc(^ — 2)c(fe; N n ) of labeled vertices in all cores of N n is in vq-n ± 
We now argue that the number total number v no of labeled vertices that lie in cores with a 
most no labeled vertices is a.a.s. in aq- pN rj^ x ^ PN I ^ N ) ^ ± Tp. This will conclude the proof. 



First note that v no = Y^k=2(k~2)c(k; N n ). Let n_ := ( TjT^STosn l ' ^ e obtain a lower 



( V(a+1) 

^£j(n)logn y 

bound on v no by keeping the summands up to n_ in the above sum: v no > ^22=2^~ 2)c(k; N n ). 
Lemma 17.11 yields that for all 3 < k < ra_ we have c(k; N n ) > (1 — e)pkaq-n with probability 
at least 1 — X^fc=2 e -0 ^ Pfen ) = 1 — So, for large n 

v no > (1 - £ )a r n J> - 2)p fc = (1 - £ )a r n J> - 2) ^ \H^o) 
(7.10) ^ ^2 r( PAr ,iVo) 

. N PnT x (pn,N ) en 

> 1 — e)arn ; ; , 

" l ' T( PN ,N ) 2' 

where the last inequality follows from the fact that the sum ^^=2 ^Pn [ xk ]T(pN, -Wo) is con- 
vergent. Using the same reasoning, we also deduce that a.a.s. 



Y,(k - 2)c(k; N n ) < (1 + e)a T n - 2) 



p k N - 2 [x k - 2 ]T(x,N ) 



(7 .11) k=2 k=2 T(PN,N ) 

- (1+£) ° T " +£n - 

Thus, to deduce a matching upper bound for -y no it is sufficient to show that a.a.s 



no 



(7.12) J2(k-2)c(k;N n ) <en. 



k=ri- 



/ \ 1/a 

Let us set n + := I , J 1 , I . Moreover, set £ = (oj{n) a+l logra) 1 '" and note that the 



oj(n) log n 

preconditions of Lemma 17.21 are satisfied with k = n+ and this choice for £. We infer that 
a.a.s. for all 3 < k < n + 

. . 1 \ 3 , 4^ [x^ 2 ]T(a;,Afo) -plT 2 
c(k, £k; N n ) < - • jPfc^ fc n, where p fcj5fc = ^ — — 

e=k °' 
Thus, by using (|7.5p we infer that there is a C > such that 

£) (fc - 2)c(fc; N n ) < C7^n f> - 2) ^"^'^ ' ^ = o(n). 
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This shows (|7.12p . and the proof is completed. □ 

8. Remaining Proofs 

8.1. Proofs of Section [2J 

Proof of Proposition ^. 1[ Firstly, let us observe that there is a natural projection of (TV + 
1) x X 2 onto Af which maps an ordered pair in (Af + 1) x X 2 to its first element, which we 
call the underlying network. Note also that for each network in Af n -2 its preimage under this 
map consists of n(n — 1) elements of (Af + 1) x X 2 . Therefore, since V is a property that 
is closed under automorphisms the proportion of networks in Af n -2 that have V equals the 
proportion of elements of (Af + 1) x X 2 whose underlying network has V. In other words, the 
probability that N n _2 £ V equals the probability that if we choose uniformly at random an 
element of (Af + 1) x X 2 then the underlying network has V. We denote by P the subset of 
(Af +1) x X 2 which is the preimage of the networks in Af n -2 which have the property V. So 
the proportion of P in (Af + 1) x X 2 is also at least 1 — f(n — 2). 

The latter probability space is mapped bijectively into (1 + e) x B n . So, in particular, P is 
mapped bijectively into a set P' C (1 + e) x B n . Let us consider the subspace e x B n which 
contains precisely half of the elements of (1 + e) x B n . We have 

1 > 1 - 2f(n - 2). 



|e x B n \ 

On the other hand, the subset e x B n is mapped bijectively into the set B+ '■= {(B,e) : 
B € B n , e € E(B)}. Therefore, if P" is the image of P' under this isomorphism, then 
Jg>l-2/ (n -2>. 

But now we are able to relate probabilities in the uniform space B£ with probabilities in 
the space B n . Note that there is a natural projection of onto B n where a pair (B,e) 6 B% 
is mapped to B G £>„,; we denote this by it. So |£>+| < nn\B n \, as every graph in B n 
has at most ktl edges. Furthermore, let P be the complement of P" in B^. Since every 
graph in B n contains at least n edges, we have |vr(P )| < \P \/n. Therefore we arrive at 

the relation < ^t^tt- Since the latter ratio is at most 2f(n — 2), it turns out that 

\On\ — |B+| J K n 

\<Q\<2Kf(n-2). □ 
8.2. Proofs of Section |H 

Proof of Theorem \4-3[ Let s = [k/m\ and abbreviate F = (1 — f /r(x, y)) 1//m . The analytic- 
ity of <?, h implies that G can be represented as 

G(x,y, f)=a Q + a m F m + ■■■ + a sm F sm + a k F k + ..., 

where the aj = ai(x,y) are auxiliary analytic functions. In particular, we have that 

(8.1) a = g(x,y,r(x ,y )), a m = -g f (x,y,r(x ,y ))r(x ,y ), and a k (x , y ) / 0. 

The definition of F implies that / = r(l — F rn ), where r = r(x, y). The (unknown) function / 
hence satisfies the equation 

(8.2) r -a = (a m + r)F m + a 2m F 2m + ■■■ + a sm F sm + a k F k + .... 
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Note that oq(xo, Do ) = 9{ x o-,Vo-, r ( x o, Ho)) = r ( x o, Ho)-, which implies that the the left-hand side 
of the above equation vanishes at (x$, yo)- Moreover, our assumptions imply that r x (xo, yo) — 
(ao)x{xo,yo) = r x {x ,y ) - g x {x ,yo,r(x ,y )) ^ 0. By applying the Preparation Theorem 
by Weierstrass we thus infer the existence of analytic functions H {x, y) and p{y) such that 
H(x ,y ) / 0, p(y ) = x and locally around (x ,y ) 

r - a = H(x,y)(x - p{y)). 

Set X = (1 — x / p{y)) l / m . Equation (18.21) is then around (xQ,yo) equivalent to 

(-H(x, y)p{y))X m = F m ((a m + r) + a 2m F m + ■■■ + a sm F^ m + a k F k ~ m + . . 



Recall that for any \x\ < 1 and a £ C we have that (1 + x) a = ^2 k>0 i^ k )x k - Using this, the 
above is equivalent to 

(-H(x, y)p{y)) l ' m X = F ({a m + r) 1 '™ + a 2m F m + ■■■ + a sm F^ m + a k F k ~ m + ... 



where the c^'s are functions given in terms of the a^'s, and in particular — m ^ a +r )(m-i)/m - 
As H(x ,y )p(y ) ^ and 

a m (x ,y ) + r(x ,y ) ^ (1 - g f (x , y , f ))r(x Q , y ) + 0, 

the above relation between X and F is locally invertible around (xQ,yo). By indeterminate 
coefficients we obtain 

k — m + 1 



/ -H(x,y)p(y) \ 1/m +1 (-H{x,y)p(y)) " +1 k _ m+l 

b = n ~J~r + tm + lX + ak ; , fc-m+2 X 

\ «m -r ' / (a m + r) m 



where the 6j's, i € {jm + 1 : 1 < j < s} are analytic functions of the c^'s. By taking the mth 
power of both sides of the above equation we obtain 

/ -H(x,y)p(y) „ {-H{x,y)p{y)) k l m k 

-Jv -)-... — CL k ui — -+- . . . . 



r(x,y) a m + r (a m + r)( fc+1 )A 

This completes the proof of (|4.4p , as it readily follows from our assumptions that the coefficient 
of X k above is ^ in a neighborhood of (xo,yo)- D 



8.3. Proofs of Sections [5] and [6l 

Proof of Lemma \5.4\ To see (|5.ip , we define a mapping from the set of vertices of the resulting 
network into the set of calls of the samplers TN, TS, TP and TH. More specifically, we map 
a vertex to the call of the subroutine where this vertex appeared for the first time. Note 
that, in fact, the image of this map is the union of the calls of TS and TH only, as these are 
the only subroutines where vertices are created. The preimage of each call of TS consists of 
only one vertex, whereas the preimage of each call of TH consists of as many vertices as the 
number of vertices of the sample from T, which is used in the call of TH. Note that once a 
vertex has been created it can never be identified with another vertex but only with a pole. 
Thus (HH]) follows. 

Similarly, (|5.2p follows from a similar mapping of the set of edges of the resulting network 
to the set of calls of the samplers TN, TS, TP and TH, where an edge is mapped to that call 
where it was created. In this case, the image of the mapping is contained into the union of 
the sets of calls of TN, TS and TP. This is the case, since in any call of TH each edge of the 
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sample from T is replaced by a network which is the result of a call of TN. and it is the calls 
of TN which may yield an edge. Also, among the calls of TS those which create an edge are 
those whose left part consists of an edge. Finally, the calls of TP which create an edge are 
those which result into a parallel network of the first type, that is, a parallel network which 
consists of an edge and a set of at least one S or Tt network. Thus, (|5.2p follows. 

Equation (|5.3p follows by mapping the set of calls of TN (apart from the initial call) to the 
union of the sets of calls of TS or TH , where a call of TN is mapped to the subroutine which 
started it. (Note that, apart from the initial call of TN, it is only these samplers that are 
able to call TN.) In particular, a call of TH makes precisely one call of TT and, conversely, 
each call of TT is made by a call of TH. Furthermore, for each edge of the sampled network 
of type T precisely one call of TN is made. This concludes the proof of (|5.3p . 

As far as fj5.4j) is concerned this is simply counting the number of calls of TP which yield a 
parallel network of the first type (if j = 1) or a parallel network of the second type (if j = 2). 

Equation (|5.5p follows with a similar argument by mapping each call of TP to the call 
of TN or TS which created it, as these are the only types among our Boltzmann samplers 
which can call directly TP. Thus (|5.5p holds and the same kind of argument works also for 
()5.7p as well as for (|5,8p . For the latter, we need the observation that there is a one-to-one 
correspondence between the calls of TH and the calls of TT. □ 



Proof of Lemma \6.1\ We start with the case j = 1. For a Po>i(^) variable X and £ > 1 we 
have E (£ x ) = n^Ej • Using Markov's inequality, this implies for any £ > 1, < r' < r, 
and e > 

(8.3) P (Y T , > (1 + e)r'E (X)) < e r ' m , where /(£) = log - (1 + e)E {X) log £. 

Note that 

/(6 = ^3T 1 and / ^ = -pr3T)i + 12 • 

The function uzfv is strictly monotone decreasing for x > 1. Together with the triangle 
inequality this implies for £ > 1 

^ j^r-T? + + e)E {x) - (E {x) + 1 + e)E (x) ■ 

Write £ = 1 + 5, where 5 > 0. We obtain by applying Taylor's Theorem 

/(£) < /(l) + f'(X)S + (E (X) + 1 + e)E (X) J 2 . 
As /(l) = and f'(l) = —eK (X), by setting 5 = 2 (i+ e +E(x)) an< ^ usm S (|8.3p we infer that 



(y r ; > (l + e)r'E(X)) < e" 



4(l+e+E(X)) 



E(y r ,) 



We can deal completely analogously with the lower tail of the distribution of Y r i - the straight- 
forward details are omitted. 

Using the above bounds we obtain readily for any < r' < r that there is a constant 
C = C(fi) > such that 



1 - e-P 



> er I < e- c ' £2r . 
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The claim for J = 1 then follows with plenty of room to spare by applying the union bound 
and using that e > . As the calculations are very similar for the case j = 2 we leave them 
to the reader. □ 
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